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In this work we present the results of a study of spectral statistics for a classically 
integrable system, namely the rectangle billiard. We show that the spectral statistics are 
indeed Poissonian in the semiclassical limit for almost all such systems, the exceptions being 
the atypical rectangles with rational squared ratio of its sides, and of course the energy 
ranges larger than I/ m ax = Ji/To, where To is the period of the shortest periodic orbit of the 
system, however L ma x — > oo when E — » oo. 

§1. Introduction 

The quantal energy spectra of classically integrable and classically chaotic sys- 
tems show a remarkable difference in their statistical properties. The former appear 
to be random, satisfying Poissonian statistics, while the latter behave as eigenvalues 
of random orthogonal (Gaussian orthogonal ensembles or GOE) or Hermitian matri- 
ces (Gaussian unitary ensembles, GUE), depending on the properties pi the system 
with respect to antiunitary transformations like the time reversal EP'cP'B 1 . 

These observations have, however, not yet been proven. There exists numerical 
evidence supporting the Poissonian model for the integrable case&. However, there 
has also been some evidence jrf the discrepancy of the Poissonian model for the case 
of the rectangle billiard© 'Q)' IIP. 

In this work we study the spectral statistics of a few integrable systems, namely 
the rectangle, torus and circle billiards. Here we present only the results for the 
rectangle billiard. For results on the torus and circle billiards and a more in-depth 
discussion of the rectangle billiard see reference 

§2. Definitions 

We deal with an ordered spectrum {Ex, E2, Ei, ...} of a given system. In 
order to compare spectral statistical properties of different systems it is convenient 
to eliminate the system dependent average density of states. This is done by a 
nonuniform mapping of the energy scale called the unfolding of spectra. If we know 
the average number of states up to a given energy N(E), which can be obtained in 
the first approximation by the Thomas-Fermi rule, then the mapping 

Xi = N(Ei) (2-1) 
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gives the unfolded spectrum {x{\ of the system whose mean level spacing is equal to 
1. 

The first measure of statistical spec- 
tral properties is the nearest neighbour 
level spacing distribution P(S), which 
gives the distribution of the distances 
between the consecutive energy levels. 
For Poissonian systems it is (after un- 
folding) 

ftwon(S) = exp(-5). (2-2) 

The next measure of the spectral 
properties is the A statistics. It is de- 
fined as 

a+L 

[N(x)-Ax-B]dx) a , 
(2-3) 

where N(x) is the spectral staircase 
function and counts the number levels 
up to energy E. This measure shows 
the average deviations of the spectral 
staircase function from the best fitting 
straight line over the interval L, and is 
therefore sometimes referred to as spec- 
tral rigidity. For systems with degen- 
eracies of zero measure, for small L it is 
always 



A(L) = — (min 
L A,B 




A(L) = L/15. 



(2-4) 



For Poissonian spectra, however, this is 
true for all L. 

Another spectral measure is impor- 
tant since all other measures can be de- 
rived from it. It is the E{k,L) distribu- 
tions, which give the probability of find- 
ing exactly k levels in the interval of the 
length L. Due to the unfolding of spec- 
tra, their maximum is close to L ~ k. 
For Poissonian systems, 

L k 

£poisson(&, L) = — exp(-L). (2-5) 
k\ 




Fig. 1. The level spacing distribution for dif- 
ferent energy windows, top to bottom: 10 6 
levels above unfolded energy 10 9 (energy 
window Oi), 10 7 levels above energy 10 10 
(O2) and 10 s levels above energy 10 11 (O3). 
In the insert the appropriate magnifica- 
tions for small level spacings are shown. 
The GOE/GUE expressions for the E(k,L) statistics can be found in the form of 
tables in the reference Slfor k = — 7, while for larger k the asymptotic Gaussian 
formulae are applicableBP . 
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§3. Results 



The billiard we studied was the rectangle. Its spectrum is in appropriate units 
given as 



2 2 

m + an , 



(3-1) 



where m, n are positive integers and a is the squared ratio of sides. 
If a, is rational, 



V 

a = -, 

q 



(3-2) 



where p and q are positive integers, the 
energy levels become integer multiples 
of the quantity If 



X 



TT 



p,q 



AJpq 



(3-3) 



in the units of the mean level spacing. 
This means that the level spacing distri- 
bution becomes a sum of delta functions 
separated by these steps. In reference 
H) we show that even the amplitudes for 
the delta peaks do not follow the Poisso- 
nian prediction and that the degeneracy 
of levels increases with energy. Connors 
and KeatingliiP proved that indeed the 
level spacing distribution for the case of 
the square a = 1 billiard goes towards 



P a=1 (S)=5(S), 



(3-4) 



with E — > oo. 

Rational a are, however, not typi- 
cal. For our studies we chose the irra- 
tional a 

suits with the value of a 



^, but also verified the re- 



V5-1 



which 



is the golden mean and which gave qual- 
itatively the same results. 

In the figure 1 we plot the level spac- 
ing distribution for different high lying 
energy windows. We see that the Pois- 
sonian statistics are followed excellently, 
with the statistical deviations remaining 
within the expected values. 

If, however, we try to include all the 
levels up to a given energy E in the plot, 





Fig. 2. P{S) for small level spacings for all lev- 
els up to energy (top to bottom) 10 7 , 10 s 
and 10 9 , with the bin size chosen such that 
the expected number of events per bin is 
about 1000. 
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Fig. 3. The A statistics for the same energy windows as in the figure 1. 

and shrink the bin sizes so that at all energies there is the same expected number 
of events in the bin, we can see in the figure 2 that there exist fluctuations around 
the Poissonian value which are larger than statistically expected and which do not 
seem to decrease with energy. We show in reference ||) that this happens due to 
the sufficiently influential rational a close to the chosen irrational value. If the bin 
size is kept constant, however, the distribution fluctuations tend to the statistically 
expected ones. 

The long range A statistics also behaves as expected, as can be seen in the 
figure ||. For L < L max , where L max = K/Tq with To being the period of the shortest 
periodic orbit, A follows the Poissonian prediction. BerryO predicted through the 
use of the periodic orbit theory that at larger L the A should saturate, as can also be 
confirmed in our plot. We may notice, however, that L max oc \[E — > oo as E — > oo. 

The E(k, L) measures for different k are plotted in figure 4. We can see that the 
Poissonian prediction is nicely followed up to k being a fraction of the L mSuX , where 
the numerical functions become narrower than the Poissonian ones. This can be 
understood by noting through the A plots that the spectrum becomes more rigid at 
higher L than expected for the Poissonian one, so we should expect exactly k levels 
to be found in a smaller ranges of L around the average value (L) = k. But, as L max 
goes to infinity with increasing energy, the Poissonian E(k, L) statistics are followed 
for all k and L in our system in the limit E — > oo. 

§4. Conclusion 

The Poissonian statistics are indeed nicely satisfied in the systems we studied, 
apart from the atypical cases of the rational squared ratio of sides of the rectangle, 
and the energy ranges greater than -L max . We show the same to be true also for the 
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Fig. 4. The E(k,L) statistics for the energy window 10 8 levels above unfolded energy 10 11 and 
k = 0..50 (top left), k = 10 3 (top right), k = 10 4 (bottom left) and k = 10 5 (bottom right) 



cases of the torus and circle billiards in reference ||). 

These studies are not important to understand the properties of the integrable 
systems only. In the case of the most general, mixed type systems, the total spectrum 
is in the semiclassical limit composed of independent contributions of the regular and 
chaotic components, with the regular component contribution having Poissonian and 
the chaotic component contribution the GOE(GUE) statistics (the Berry-RobniklllP 
theory), so a good understanding of the individual components is necessary in order 
to understand the whole picture. For a study of mixed type systems see reference 



14 ) and the short version in reference [I 
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